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ABSTRACT 

We study via numerical N-body/SPH chemistry simulations the effects of primordial non- 
Gaussianities on the formation of the first stars and galaxies, and investigate the impact of 
supernova feedback in cosmologies with different /nl- Density distributions are biased to 
higher values, so star formation and the consequent feedback processes take place earlier in 
high-/NL models and later in 1ow-/nl ones. Mechanical feedback is responsible for shock- 
ing and evacuating the gas from star forming sites earlier in the highly non-Gaussian cases, 
because of the larger bias at high densities. Chemical feedback translates into high-redshift 
metal filling factors that are larger by some orders of magnitude for larger /nl, but that con- 
verge within one Gyr, for both population III and population II-I stellar regimes. The efficient 
enrichment process, though, leads to metallicities > 10~^ by redshift ~ 9, almost inde- 
pendently from /nl- The impact of non-Gaussianities on the formation of dark-matter haloes 
at high redshift is directly reflected in the properties of the gas in these haloes, as models 
with larger /nl show more concentrated gas profiles at early times. Non-Gaussian signatures 
in the gas behaviour are lost after the first feedback takes place and introduces a significant 
degree of turbulence and chaotic motions. Despite this, our results support the idea that non- 
Gaussianities can be imprinted in the gaseous and stellar features of primordial structures in 
the high-redshift Universe. 
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> 1 INTRODUCTION 

Among the sever al outstanding problems of modern Astrophysics 
° and Cosmologv (iBarkana & Loebl 1200 ll : ICiardi & Feiwal l2005l : 
5^ \ iBromm & Yoshidal 12011 ^, the status of the early Universe is cer- 
tainly a very debated one. The formation of primordial struc- 
tures is supposed to be strongly influenced by the cosmologi- 
cal initial conditions they originated from, and it is widely ac- 
cepted that all the visi ble objects derive from the growth of mat- 
ter perturbations (e.g. "Ounn & Gott* 'l972'; 'White & Rees ' l978|; 
Peebles. P. J. E.ll993: Peacock. J. A. 1999; Coles, P & Lucchin, R 



200HPeebles & Ratral l2003). that developed shortly after the Big 



Bang, during the era of inflation. These perturbations have grown 
in time in an expanding Universe and they have assembled into the 
galaxies, galaxy groups, and galaxy clusters observed today. 
Acc ording to recent determ inations of the cosmological parameters 
(e.g. iKomatsu et al.ll201ll) . the Universe is composed by ^ 30% 
of matter and for the remaining ~ 70% of an unknown term at- 
tributed to th e so-called cosmo logical constant. A, or dark energy 
(see also e.g. Mai o et al.ll200 ^, for N-body/SPH chemistry simula- 
tions in dark-energy cosmologies and the effects on baryonic struc- 
ture evolution). More precisely, the present matter contributions 
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to the cosmic density are (IKomatsu etalll201l h ^^o,m 0.272, 
^^o,A = 0.728, Qo,b = 0.044, for matter, cosmological constant, 
and baryons, respectively. The cosmic equation of state parameter 
is consistent with w — —1, the observed spectral index of primor- 
dial fluctuations is n = 0.96, and the power spectrum normaliza- 
tion is given by a mass variance within 8 Mpc//i-sphere as = 0.8. 
Structure formation depends strongly on the initial density 
fluctuations imprinted on the primordial m atter distribution 
dPress & Schechter 1974l : lsheth & Tormen|[T999) . The standard as- 
sumption on the distribution of density fluctuations in the Uni- 
verse is based on an episode of dramatic size growth of the uni- 
verse, roughly 10~^^ s after the Big Bang, during which the 
seeds of the prese nt-day structures formed dStarobinskvl I198Q| : 
lOuthl Il98ll : iLindd [1990). These models of inflation predict in 
general that the overdensity ^ is a Gaussian random variable 
with variance fully determined by the underly i ng power- spectrun i 
(e.g. [K omatsu 2010; Casapo nsa et al.l l201ll : ICurto et al.l l201ll : 
iBruni et al.l 120111 . and references therein). The general consen- 
sus on the Gaussianity derives mainly from the central limit 
theorem. However, o bservational determinations (see Table in 
iMaio & Iannuzzi|[201 1) show evidence for deviations from Gaus- 
sianities that leave room for non-Gaussian investigations, as 
wen (P eebles 1983; Viel et al. 2009; Desiacaues & Seliak 2oTqI: 
ILo Verde & Smith i2011i : iDesiacques et al.. .2011. : .D'Amico et al.l 
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I2OIII : iHamaus etal1l201lh . The effects of non-Gaussianities are 
expected to play a role mostly for the high-sigma density fluctua- 
tion s (e.g. Grinstein & Wise 1986: Kovama et al.'l999^: lzaldarriagal 
I2OOQ : Wagner et al. 2010; Lo Verde & Smith 2011), and, thus very 
early structures should result somehow affected by them jMaiol 
l201lh . In particular, due to the sensitivity of the gas cooling ca- 
pabilities to the underlying matter densities, the initially skewned 
non-Gaussian features could be reflected by the earlier collapse 
of molecular gas, and theoretically influence the formation epoch 
of the first stars and galaxies, as pointed out via numeric al N- 
body/SPH chemistry simulations bv lMaio & lannuzzil (l201lh . Be- 
yond the formation redshift, also the consequent feedback mecha- 
nisms could play an important role in ejecting material from the star 
forming regions or in polluting the surrounding environment. Be- 
cause of the lack of relevant studies dealing with non-Gaussianities 
and feedback mechanisms, these are still open questions which we 
will discuss and investigate throughout this work. 
In particular, it is possible to distinguish into mechanical, chemical, 
and radiative feedback (for an extensive review of the many pos- 
sible mechanisms, see e.s. ICiardi & Ferr^l2005L and references 
therein). The first class includes all the different phenomena re- 
lated to mass or energy deposition into the cosmic gas from star 
formation and stellar evolution (i.e. shocks, blowout, blow-away, 
etc.); the second one comprises essentially the effects of chemical 
changes in the gas composition (i.e. metal enrichment and conse- 
quent stellar population transitions); and the third one covers the 
aspects linked to radiation emitted by cosmic sources (i.e. pho- 
toionization or photodissociation of molecules and atoms, gas pho- 
toheating, cosmic reionization, etc.). 

We will mainly consider mechanical and chemical feedback from 
first structures, both from population III (popIII) and from popula- 
tion II-I (popll) stars. The transition between these two regimes 
is determined by a critical metallicity, Zcrit which, according 
to different authors dSchneider et aDl200l l2006l : E romm &Loebl 
l2003h . is estimated to be around ~ 10"^ - lO'^Zc;). It has also 
been previously invest igated in details (e.g. iTornatore et al]|2007l : 
lMaioetal.ll20ldl201lh with the help of numerical simulations fol- 
lowing the chemical evolution of the cosmic gas and metal pollu- 
tion from stars. The substantial distinction between the popIII and 
the popll regime is the stellar initial mass function. While in the 
latter case, it is known to be Salpeter-like, in the former case it is 
unknown, even if expected to be top-heavy, for the incapability of 
pristine gas to cool down to very low temperatures, in presence of a 
high CMB floor (a few hundreds K at ^ :^ 10). However, fragmen- 
tation of primordial gas clouds and formation of popIII star with 
masses below ^ 10^ can still be possible (as shown by e.g . 
Yoshidal 12 006": 'Yoshid a et all l2007l : ICampbell & Lattanzid l2008l : 



Suda & F uiimoto 2010). The impacts of different assumptions on 



the primordial IMF, yields, and supernova ranges have already been 
largely studied by Maio et al. (2010), thus we will not go into the 
details here and will simply assume a top-heavy IMF. Stars with 
masses larger than ^ 100 M© are short-lived (up to ^ 10^ yr) 
and the ones in the range [160, 240] M© die as pair-instability SN 
(PISN) and are able to pollute the surrounding medium by ejecting 
large amounts of metals. 

In the present work, we will focus on the impacts of mechanical and 
chemical feedback in the primordial Universe and address their role 
in presence of non-Gaussian initial conditions. The paper is struc- 
tured as follows: after presenting in Sect.[2]the simulations used, in 
Sect. [3] we will discuss the main results related to the cosmologi- 
cal effects of mechanical feedback (Sect. 13.11 ) and chemical feed- 
back (Sect. l3.2l ), and we will also show the consequences for single 



haloes over cosmic time fSect. l3.3l ) and the implications of primor- 
dial streaming motions (Sect. [34l ). In Sect. |4] we will summarize 
our findings and conclude. 



2 SIMULATIONS 

The simulations considered here were firstly described in 
iMaio & lannuzzil (I2OIII) . who performed a large set of runs with 
different box sizes and resolutions. Since we want to study in detail 
the joint non-linear effects of feedback mechanisms and primor- 
dial non-Gaussianities, we will focus on the high-resolution simu- 
lations, having a box size of 0.5 Mpc//i, and an initial gas parti- 
cle mass of ^ 40 M© / h. Local non-Gaussianities are included by 
adding secon d-order perturbations to the Bardeen gauge-invariant 
potential (e.g. lSalopek & Bondl[T99Q i): 



^ = $L + /nl < >] 



(1) 



with $L the linear Gaussian part, and /NLthe dimensionless cou- 
pling constant controlling the magnitude of the deviations from 
Gaussianity. Observational constraints on /nl suggest values be- 
tween ~ — 100 (a complete table of observational determinations 
is given in lMaio & Iannuzzill20 1lV so we will focus on the cases 
/nl=0, /nl=100, and also on /nl=1000 for sake of comparison. 
The simulations were performed by u sing a rnodified version of 
the parallel tree/SPH Gadget-2 code (ISpringel| l2005h . which in- 
cluded gravity and hydrodynamics, with radiative gas c ooling both 
from molecules and atomic transitions (according to iMaio et aD 
l2007h . multi-phase model ( Springel & Hernguist 2003h for s tar for- 
mation, UV background radiation jHaardt & Madaul 199^ . wind 
feedback ("Sprin gel & HernQuistll2b03l : lAguirre et al.ll200lh . chemi- 
cal network for e~,H, H+,H~,He, He+, He++, H2, H+, D, D+ 
HP, HeH+ (e.g.lYoshida et al. 12 003: Maio et al. 200d l2007ll2009l: 
lMaioll2009l :l Maio et al .1120 id and references therein), and metal 
(C, O, Mg, S, Si, Fe) pollution from popIII and/or popll stellar 
generatio ns, ruled by a c ritical metallicity threshold of Zcrit = 
10""^ Z© (ITornatore etal.l f2007: Maio et al. 2010, 2011). The cos- 
mological parameters are fixed by assuming a concordance ACDM 
model with matter density parameter r2o,m = 0.3, cosmological 
density parameter Q^q^a — 0.7, baryon density parameter Qo,b = 
0.04, expansion rate at the present of Ho = 70 km/s/Mpc, power 
spectrum normalization via mass variance within 8 Mpc//i radius 
sphere as = 0.9, and spectral index n — 1. We consider a top- 
heavy stellar initial mass function (IMF) with mass range [100, 



500] M0 for the population III regime (different c ases are dis - 
cussed in 'Maio et al.(l2oTQl : iMaio & lannuzzil 1201 ll : iMaiol 1201 ih . 
and a Salpeter IMF with mass range [0.1, 100] M© for the pop- 
ulation II-I regime. 

A friend-of-friend (FoF) algorithm with comoving linking length 
of 20 per cent the mean inter-particle separation, is applied at post- 
processing time to find the formed cosmic objects, with their dark, 
gaseous, and stellar components. 



3 RESULTS 

In the following we will present results related to the inter- 
play of mechanical feedback and chemical feedback with non- 
Gaussianities (Sect. 13.11 and Sect. 13. 2K and the implications for 
early cosmic structures (Sect. [33] and Sect. 13.41 ). 
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3.1 Mechanical feedback and gaseous properties 



We begin our discussion by commenting on the thermodynamical 
properties of the cosmic gas and the feedback mechanisms related 
to stellar evolution and supernova explosions. 



3.1.1 Distributions 

As already mentioned, non-Gaussianities play a role on the primor- 
dial matter distribution. Therefore, we start our investigations by 
studying the gas mass distribution for the different /nl scenarios. 
We note that star formation sets in when the gas reaches a density 
of - 10^ cm . At this point, feedback mechanisms become ac- 
tive, too. In Fig.[TJ we plot the cumulative mass fraction of gas hav- 
ing density larger than a given value, at redshift z ^ 23, z :^ 13, 
and z ^ 9 for /nl=0, /nl=100, and /nl=1000. At z - 23 (left 
panel), most of the gas is still in low density environment, with only 
a small fraction of — 10 ~^ undergoing first collapsing episodes, at 
number densities n > 0.1 — 1 /i^cm~^. The /nl=0 and /nl=100 
models are almost identical, while the /nl=1000 case shows a 
larger tail at the high-density end, demonstrating that H2 and HD 
molecules have been more efficiently formed and the molecular 
gas content has already reached fractions of> 10~^ — 10~^. 
This allows the gas to cool and condense more rapidly and to 
reach n 10^ /i^cm~^, while in the other models densities of 
n — 1 h^cm~^ are barely achieved. At later times, when ^ :^ 13 
(central panel), the density range is equally covered for all the /nl, 
and the contribution from clumped regions increases of a few or- 
ders of magnitude. The /nl=1000 cosmology preserves some sig- 
natures of the primordial distribution, and these are reflected in a 
factor of 2 — 3 in the higher distribution for n > 0.1 h^cm~^ . 
Finally, at z ^ 9 (right panel), the behaviours converge with some 
per cent of the gas collapsing into the dense star forming regions 
and almost erasing any evidence of non-Gaussianity. Residual con- 



tributions persist around n 10 
though. 



10^ h^cui-^ for /nl=1000. 



3.1.2 Global evolution 

In order to discuss more in depth mechanical feedback and its im- 
plications in Gaussian and non-Gaussian cosmologies, we study the 
global behaviours of gas density over cosmic time. 
In Fig. [21 we plot the redshift evolution of gas densities, by focus- 
ing on the average density, the minimum density, and the maximum 
density reached in the runs. 

We note that the mean mass density, i.e. the ratio between the to- 
tal mass in the box and the simulated volume does not give any 
significant information on the different models, since it evolves 
oc (1 -\- independently from /nl- Thus, in order to see local 
signatures of non-Gaussianities on the gas behaviour and on its 
collapsing phases, we consider the average value of the densities 
of each gas particle, computed in the SPH fashion by smoothing 
the surrounding mass over the SPH kernel (in different words, the 
average we are showing is the average of the density PDF in each 
of the simulations). In this way we can easily identify effects from 
the underlying /nl on structure growth. In fact, larger /nl will en- 
hance structure formation and hence hydrodynamical densities will 
increase more rapidly. 

According to the figure, the average densities follow the cosmic av- 
erage evolution fairly well and decrease smoothly down to redshift 
z ^ 20 (corresponding to — 10 ~^ h^cm~^), when gas collapses 
and star formation episodes become effective, and boost the mean 
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Figure 2. Redshift evolution of the gas densities in the simulations for 
/nl=0 (solid lines), /nl=100 (dotted lines), and /nl=1000 (dashed lines). 
The data related to maximum, average and minimum values for each model 
are plotted from top to bottom, as indicated by the labels. 



values by 2 orders of magnitude. This is well seen in the devi- 
ation from the decreasing trends from high to low redshift which 
signs the departure from the mean cosmic value and the onset of 
star formation. In fact, at z < 20, when first highly non-linear 
structures form, densities increase locally of several orders of mag- 
nitude, and dominate over the typical values of the other particles. 
Basically, the different averages for different /nl at lower redshifts 
reflect the different non-linear behaviours of the very first collaps- 
ing objects. 

The universes with /nl=0 and /nl=100 have very similar aver- 
ages, whereas the model with /nl=1000 shows earlier deviations 
because of the earlier collapse phases undergone by the gas. The 
trend of the /nl=1000 case is easily understood when looking 
at the maximum densities sampled by the simulations. Both for 
/nl=0 and /nl=100 densities at early times (z > 20) are around 
1 — 10 h^cm~^, instead, for /nl=1000 they are systematically 
higher of a factor of a few and, thus, can grow faster, due to po- 
tentially enhanced cooling instabilities. This is well visible during 
the first collapse, around z ^ 20, when the maximum density in- 
creases exponentially of ^ 2 orders of magnitude in ^ 20 Myr, 
and the discrepancies between /nl=1000 and /nl=100 or /nl=0 
are particularly evident. The minimum densities do not show sig- 
nificant modifications with non-Gaussianities. 
In all the models, early evolution is characterized by a medium 
which is quite cold (at temperatures of a few hundreds Kelvin) 
and that is cooled predominantly by primordial molecules. Shock 
heating of gas while it is falling into the dark-matter potential 
wells causes temperature increase up to ^ 10^ — 10^ K. When 
the first star formation episodes take place (at z ^ 20) in the dens- 
est regions, stellar feedback, through PISN/SN explosions, rapidly 
brings the temperatures to ^ 10^ — 10^ K, several orders of magni- 
tude hotter than the cooling medium. Hot gas shocks the surround- 
ing material and pushes the gas to lower-density regions, as well (it 
is the simultaneous ongoing gas collapse in different sites keeping 
high the maximum densities of Fig.|2]). These sudden changes allow 
us to detect the earlier structure evolution mainly in the /nl=1000 
universe. These stages are very short, though, and in some tens of 
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Figure 1. Gas cumulative probability distributions, /(> n), as a function of number density, n, at redshift, 2; ~ 23 (left), z ~ 13 (center), and z ~ 9 (right), 
for /nl=0 (solid lines), /nl=100 (dotted lines), and /nl=1000 (dashed lines), as indicated by the labels. 



Myr the different values in three models converge and the discrep- 
ancies fade away. 

3.2 Chemical feedback and metal enrichment 

At this point, we discuss the main results related to the chemical 
feedback and its effects, in particular metal pollution at early times. 
A pictorial representation of the pollution events in the different 
boxes is given in Fig.O where we display the metallicity maps for 
/nl=0 (left), /nl=100 (center), and /nl=1000 (right), at z = 15 
(upper row) and z — 9 (lower row). These immediately show some 
differences among the various cases, mostly at 2; = 15, when the 
first spreading events are taking place. For all the cases, enrichment 
up io Z/H ^ 10~^ is detected, but it is more limited for /nl=0 
and /nl=100 than for /nl=1000. At redshift z = 9, the amount of 
metals spread is comparable, but, as we will better quantify in the 
next sections, in the /nl=1000 there is slightly more enrichment 
that results in a larger filling factor. 

3. 2. 1 Phase distributions 

In Fig. 151 we show phase diagrams at different redshift for the en- 
riched particles, color-coded according to their probability distribu- 
tion. Early star formatio n episodes, domin ated by massive, short- 
lived population III stars (iMaio et aDl201(t) quickly enrich the sur- 
rounding medium, from the higher to the lower densities. The tim- 
ing of these episodes is affected by the adopted /nl, though, mostly 
at very high redshift. Differences for /nl=0 and /nl=100 are quite 
small, but they become much more evident for /nl=1000. In this 
latter case, in fact, the densities (see also next sections) are strongly 
skewned towards higher values, thus gas cooling, condensations, 
and fragmentation is more e nhanced and lead to an earlier onse t 
of star formation (Fig. 4 in llMaio & lannuzzil [20T1I : lMaidl201lh . 
Indeed, first star formation events are already detected at 2; ^ 23 
(when the cosmic time is about 144 Myr) for the /nl=1000 case, 
while they appear only at z ^ 19 - 20 (^ 45 - 32 Myr later) 
for /nl =0-100. Given the rapid evolution of primordial massive 
stars (< 20 Myr), metal pollution in the /nl=1000 case has al- 
ready taken place, when the first stars are formed in the models with 
/nl=0 and /nl=100. By comparing the time sequence of the pol- 
lution patterns for /nl=0, /nl=100, and /nl=1000, we note that at 
2; ^ 18 they show quite different stages. For /nl=0 and /nl=100, 
the first metals are just being ejected from newly born stars in high- 
density environments (and the first heated particles are also visible 



at T ^ 10^ K), while for /nl=1000 the enrichment process is in a 
more advanced stage, with enriched materials spread much further 
out, and reaching also very underdense regions. 
However, the following spreading episodes mitigate the effects of 
non-Gaussianities and by z ^ 15 (i.e. in roughly 60 Myr) the metal 
distributions in the phase space become quite similar. 

3.2.2 Probability distributions 

The trends for metal spreading are supported by the probability dis- 
tribution functions displayed, for sake of completeness, in Fig. |5] 
For these calculations we define the overdensity variable: 

5 = p/pcr, (2) 

with p gas density, 

the expansion parameter 

if = ifo [(1 - nO,tot)(l + Zf + Qo,m(l + Zf + Qo,a] , 

(4) 

and Qo.tot ^^o.m + ^^n a (iPeebles. P. J. E.lll993l : IPeacock. J. a] 
ll999l : IColes. P & Lucchin. F.II2O O2'). 

In Fig. 121 we consider metallicity, Z, and overdensity, ^ - as defined 
in the previous eq. (IJ}, (|3]l, and Q. The plots show the joint proba- 
bility distributions P{Z,6) (left column) and the conditional prob- 
ability distributions P{Z\6) (central column) and P{6\Z) (right 
column), for /nl=0 (upper row), /nl=100 (central row), and 
/nl=1000 (lower row), at redshift z = 9.00. In the plots it is 
evident the efficient enrichment process spanning a large range of 
over-densities and metallicities. Metals are produced at high densi- 
ties (6 > 10^) and ejected by wind feedback into the surrounding 
pristine regions. The high metal yields allow most of the involved 
gas to be rapidly polluted to metallicities Z > Zcrit and to pass to 
the popII-I regime, almost independently of /nl, as well demon- 
strated by all the probability distributions. The differences surviv- 
ing at this redshift are difficult to detect. The joint probability dis- 
tributions (left column) for /nl=0 and /nl=100 are basically iden- 
tical, with most of the gas atlO~^ < Z < 10~^ and some residual 
gas below ^ 10~^. Also the /nl=1000 case is very close to the 
previous ones, and no significant differences are found. 
The conditional probability distributions (central and right 
columns) show some slight deviations of the /nl=1000 case with 
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respect to /nl=0 and /nl=100, but statistically they are not rele- 
vant, and all the cases seem to have gas enriched well above Zcrit 
and spread down io 5 ^ 10~^ — 10~^. 

Comparing to the findings regarding higher redshift (previous sec- 
tion), we can state that differences in the enrichment episodes for 
different non-Gaussian cosmologies are occurring mainly in pri- 
mordial epochs, and therefore, the resulting metal filling factor 
could have signatures depending on the various scenarios. 



3.23 Filling factors 

As very first enrichment episodes can have distinct impacts in 
Gaussian and non-Gaussian models, it is relevant to quantify the 
metal filling factor at early times. 

In Fig.[6l we show the redshift evolution of the metal filling factors, 
fv, for the three cases considered, /nl=0, /nl=100, /nl=1000. 
We define fv according to: 



fv 



V 



(5) 



with i and j integers running over the number of the enriched par- 
ticles and of all the SPH particles, respectively, mz particle metal 
mass, m particle total mass, pz metal mass density, p total-mass 
density, and V simulation volume. We perform the calculation for 
the whole simulations and, additionally, by distinguishing for the 
two stellar populations (popIII and popII-I). In Fig. [6] we plot the 
evolution of the global filling factors and of the popIII filling fac- 
tors (lower lines, denoted by fzKZ^^n )■ The redshift interval (from 
z 2?> io z 9) covers a period of ^ 400 Myr, from when 



the Universe was 140 Myr old, to when it was ^ 540 Myr 
old. The results show well the distinct behaviours of the different 
models at early times, with an enrichment process which is ear- 
lier in the /nl=1000 case and later in the /nl=0 case. At redshift 
z ^ 19, the pollution levels differ of about two orders of mag- 
nitude between /nl=0 and /nl=100, and between /nl=100 and 
/nl=1000. The continuous spreading events alleviate the delay be- 
tween /nl=0 and /nl=100 to a factor of a few by z ~ 16, and 
lead to full convergence by z ^ 12. Also the advanced status of the 
enrichment process in the /nl=1000 model is gradually reduced 
to about one order of magnitude hy z ^ 14 and achieves almost 
convergence at z ^ 9. The trends for the popIII regimes are similar 
and their contribution to the total filling factor always drops down 
to 1% between z 16 and z 9 (i.e. in 300 Myr), de- 
spite the different /nl values adopted. The transition from the pri- 
mordial popIII regime to the following popII-I regime is regulated 
by the stellar lifetimes and initial mass function of the first popu- 
lations dMaio et al.l 12010), so different assumptions on the popIII 
IMF would change the relative contributions of the popIII regime 
and the onset of metal pollution, but would not substantially mod- 
ify the conclusions on the effects on non-Gaussianities. Indeed, the 
different behaviours are consequenc es of the earlier onset of star 
formation in the higher- /nl model dMaio & Iannuzzill201ll : iMaiol 

laonl). 



3.3 Implications for cosmic haloes 

For a further investigation, we study the implications of feedback 
mechanisms on the cosmic haloes in Gaussian and non-Gaussian 
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Figure 4. Phase distributions for metal enriched particles at redshift z = 18.01 (upper row), z = 15.00 (central row), and z = 9.00 (lower row), for /nl=0 
(left column), /nl=100 (central column), /nl=1000 (right column). The x-axes are comoving densities, and the y-axes are temperatures, in logarithmic scale; 
the colors refer to the probability distributions of the enriched particles. 
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Figure 6. Metal filling factors as function of redshift in the 0.5 Mpdh 
side boxes (left), and 0.5 Mpc/h side box (right), for /nl=0 (solid lines), 
/nl=100 (dotted lines), and /nl=1000 (dashed lines). Individual contribu- 
tions from popIII regions (0 < Z < Zcrit, bottom lines) are also shown, 
as indicated by the legends). 



models. As mentioned earlier, the haloes and their respective dark, 
gaseous, and stellar components are found by use of a FoF algo- 
rithm. 



3. 3. 1 Cosmic gas within formed structures 

In Fig. [7] we plot the redshift evolution of the gas mass in the largest 
objects and the corresponding gas fractiorQ • Also in this case, for 
/nl=0 and /nl=100 there are no evident differences, while for 
/nl=1000 there is an earlier formation of the structures, consis- 
tently with the previous sections. 

The mass growth is quite different at very early epochs and in fact, 
in the /nl=0 and /nl=100 cases 10^ Mq / h objects are formed 
at z ^ 30, while in the /nl=1000 case they are formed already at 
z ^ 40, about 40 Myr earlier, and are almost one order of magni- 
tude larger at z ^ 30. Such huge discrepancy shrinks over the time 
and is reduced down to a factor of ^ 2 at z ^ 20, when the gas 
masses contained in the primordial haloes are ^ 10^ Mq/Zi, and 
completely disappears at redshift z < 15. 

We note the gas that gets accumulated within early structures has a 
direct connection to the growth of the dark-matter haloes and thus 
reflect the original differences in /nl- The baryonic build-up, with 
the catch-up of the dark-matter structures, is well visible at high 



We note that at high redshift the number of particles per halo is rather 
small, and so the error on the estimated quantities is large. 
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Figure 5. Joint probability distribution P(Z,6) (left column) and the conditional probability distributions P(Z\6) (central column) and P(S\Z) (right 
column), for /nl=0 (uppe r row), /nl ^IQQ (central row), and /nl=1000 (right row), at redshift z = 9.00. The vertical lines refer to the density threshold for 
star formation (for details iMaio & Iannuzzi.2011.) . while the horizontal line correspond to the critical metallicity, Zcrit = lO^^Z©. 
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Figure 7. Mass (top) and fraction (bottom) evolution of the gaseous compo- 
nents for the biggest halo as a function of redshift, for /nl=0 (solid lines, 
cross symbols), /nl=100 (dotted lines, rhombic symbols), and /nl= 1000 
(dashed lines, triangular symbols). 



redshift, when gas is falling into the first dark-matter haloes and 
starts cooling. The first SN bursts and feedback effects are then the 
main responsible for gas expulsion from star formation events, and 
for spreading material into the low-density environments by winds. 
These processes contribute to settle down a self-regulated regime in 
which infalling gas is gradually converted into stars and feedback 
mechanisms control the amount of material which is expelled. 



3.3.2 Halo profiles 

To conclude, in Fig. [8] we study the effects on the profiles of the 
largest haloes at different epochs. 

Broadly speaking, the temperature profiles (upper panels) are char- 
acterized by high central values (^ 10^ — 10^ K), where most of 
the star formation is taking place, and lower ones in the periph- 
eral regions, where the gas is being shock-heated to ^ 10^ — 10^ K 
and subsequently cools down via molecular emissions to some hun- 
dreds Kelvin (e.g. the drops at ^ 20 — 30 pc//i at z = 15, or at 
50 - 100 pc//i at 2; = 9). At redshift z = 15, during the first 
stages of star formation, non-Gaussianities are slightly visible for 
/nl=1000 and the only effect is a larger temperature of a factor of 
^ 2 — 3 at a distance of 1 kpc//i from the center. This is related 
to the larger halo mass, which results in higher shock heating tem- 
peratures, and in an indirect effect of /nl- At redshift z = 9, star 
formation and feedback mechanisms have been going on for almost 
half a Gyr, the turbulent state of the medium heavily dominates the 
profiles, and it becomes impossible to identify trends correlated to 
/nl. 

Density profiles (lower panels) have a very regular, smooth shape 
at z = 15, and the most clumped regions correspond to the coldest 
temperatures. At this epoch, the larger concentrations in the core 
of the structures could be associated to non-Gaussianities, mostly 
for /nl=1000. Precisely, this is a more directly linked to the ear- 
lier build-up of halo masses in the /nl=1000 model. Later on, the 
cumulative development and propagation of shocks from SN ex- 
plosions determine a complex interplay between infalli ng cold ma- 
terial and expanding hot gas (see a detailed treatment in lMaio et al.l 
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Figure 8. Mass-weighted temperature, Tmw (upper panels), and mass den- 
sity, pgas (lower panels), profiles as a function of physical radius, r, at red- 
shift z = 15 (top) and z = 9 (bottom), for /nl=0 (solid lines), /nl=100 
(dotted lines), and /nl=1000 (dashed lines). 



ll), so, the continuous gas compressions and rarefactions (as e.g. 
at radii of ^ 100 pdh) erase the signal from different /nl and wash 
it out by 2; 9. 



3.4 Feedback, non-Gaussianities and primordial bulk flows 

We conclude by briefly showing the effects of primordial streaming 
motions on the formation of primordial structures and metal enrich- 
ment from first stars, in Gaussian and non-Gaussian models. As an 
example, in Fig.|9l we plot the metal filling factor for the high res- 
olution simulations with /nl=0 re-run with initial primordial gas 
streaming motions at decoupling of ^;b=0, 30, 60, 90 km/s. Such 
primordial steaming motions determine gas bulk flows on ~ Mpc 
scales and hinder gas condensation in mini-haloes with masses 
below - 10^ Mrr^ /h (iTseliakhovich & Hiratal |2oTqI : iMaio et all 
l201lh . because larger potential wells are needed to equal the ad- 
ditional contribution to gas kinetic energy. As a result, the onset of 
star formation and consequent feedback mechanisms are delayed to 
lower redshift, when the haloes have grown enough to trap the gas 
and let it cool and condense. 
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Figure 9. Metal filling factors as function of redshift for /nl=0 and primor- 
dial gas bulk velocity v\^=0 km/s (solid lines), 'Ub=30 km/s (dotted lines), 
v\y=60 km/s (dashed lines), and V]^=90 km/s (dot-dashed lines). Individual 
contributions from popIII regions (0 < Z < Zcrit , bottom lines) are also 
shown, as indicated by the label. 



4 DISCUSSION AND CONCLUSIONS 

In this work, we have discussed the principal implications of feed- 
back effects on non-Gaussianities, by using the N-body/SPH chem- 
istry simulations, presented by Maio & lannuzzi ( 2011). 
Besides gravity and hydrodynamics dSpringell l2005h . the runs 
include radiative gas coo ling both from m olecules and atomic 
transitions (according ^ to Maio et "Zl l2007h . multi-phase model 
("Sprinsel & Hernquist 2003*) for star formation, UV background 
(Haardt & Madau 1996), wind feedback ( Springel & Hernquist 
[2003l : lAguirre et al.ll200ll) . chemical network for e~ H, H+, H~ 
He, He+, He++ H2, H^, D, D+, HD, HeH+ (e.g. ' YoshidaetaP 
I2OO3I : iMaio et all 12004 1200I and references therein), and metal 
pollution from popIII and/or popII-I stellar generat ions, ruled by a 
critical metallicity threshold of Zcnt = 1Q~^ jTornatore et al.l 
l2007l : lMaio et al.l20ldl2oTTI : lMaiol201lh . Local non-Gaussianities 
are included by means of the /nl parameter, here chosen to be 
/nl=0, /nl=100, or /nl=1000, in order to explore thermodynam- 
ical and statistical properties of gas evolution and of the impacts 
of feedback effects on cosmic structures in different non-Gaussian 
universes. 

Since non-Gaussianities are based on the primordial matter density 
field, it is natural to start our study from the gas cumulative distribu- 
tion.The different models show deviations at large densities, consis- 
tently with the biased initial conditions in larger- /nl models, and 
at high z, first-molecule formation and catastrophic gas collapse are 
strongly enhanced by larger non-Gaussian parameters. The global 
thermodynamical evolution of the gas shows accordingly signifi- 
cant deviations from the Gaussian case mostly at high densities, for 
larger /nl- As a consequence, star formation and heating due to 
PISN/SN explosions take place earlier and chemical feedback, as 
well. Metal spreading preserve some information on the original 
/nl values with wider areas involved by the pollution process for 
/nl=1000 and smaller areas for /nl=100 and /nl=0. The differ- 
ences are more visible during primordial epochs (z > 15) and less 
and less at later times (z < 9), as well demonstrated by phase and 
probability distributions. The enrichment process is very efficient in 



polluting the medium and even underdense regions (6 < 1), achiev- 
ing Z > Zcrit in few tens of Myr. The metal filling factor initially 
(at z > 15) shows discrepancies up to ^ 2 orders of magnitude, 
but successively the trends rapidly converge within a factor of a few 
(at redshift z ^ 9), for both popIII and popII-I regions. By looking 
at the largest-halo mass evolution, it emerges that not only early 
structures form earlier for larger /nl, but also that their potential 
well growth can trap more gas, because of the larger halo masses. 
At lower redshift, the trends converge, regulated by star formation 
processes. The most-massive halo temperature and density profiles 
confirm that and suggest the fundamental reason why signatures of 
non-Gaussianities in the luminous structures are washed out: i.e., 
mechanical and chemical feedback strongly perturb the surround- 
ing medium, the gas behaviour becomes soon dominated by shocks 
and turbulence and, thus, looses memory of its original distribution. 
We stress that, despite the small box size (which prevents extremely 
large objects to be formed), the high resolution 40 Mq/H per 
gas particle) of the simulations enabled us to check for the first time 
the effects of non-Gaussianities on quite tiny scales (up to a few 
parsecs), at high redshift. In general, there is a delay of ^ 50 Myr 
in /nl=0 and /nl=100, with respect to /nl=1000. These effects 
are visible at very early times, without the need for implement- 
ing further non-Gaussian corrections (like the on es related to .g'NL, 
or tnl parameter s, as also discussed in detail bv iMaio & lannuzzi 
l201ll : lMaioll20TTI) . 

Whereas, the inclusion of non-linear effects that imply primor- 
dial gas streaming motions with velocities Vh ^ 30 km/s (e.g. 
iTseliakhovich & Hiratall2010l : lMaio et aLbOllD could sli ghtly alter 
the in itial phases of star formation (as demonstrated bv ^Maio et alj 
201 ih and add some degeneracies with the /nl effects (see 



Maio & lannuzzil '201 1). Substantial differences on the feedback 



impacts are not expected, though. 

Some changes would appear in case of a different popIII IMF, 
which, at the present, is basically unknown, and would be mainly 
related to the diff erent lifetimes of the different stellar masses sam- 
pled (see also Maio & lannuzzil I20 lib . For the top-heavy popIII 
IMF considered here, the mass range was between IOOM0 and 
5OOM0, and the first stars to die and pollute the medium were the 
PISN, in the range [140, 260] M©, after at most 2 Myr. For low- 
mass popIII IMF the first stars exploding would be ^ lOM© SNe, 
after ~ 10 — 100 Myr. This means that the stellar evolution process 
of low-mass (< IOOM0) stars would slow down the activation of 
feedback mechanisms and postpone the spreading events of some 
tens or a hundred Myr. 

In general, the main differences among the various cosmologies are 
detected at high redshift, when luminous structures start their life. 
This happens because larger /nl parameters determine a density 
field which is biased to higher densities, lead to larger most mas- 
sive haloes, and, as a consequence, imply shorter cooling times for 
the amount of gas that accumulates in the center. Since stars form 
earlier, metal enrichment takes place earlier, too. 
Thus, all these processes can be seen as derived effects of the non- 
Gaussian initial conditions. In particular, the essential point is the 
the formation time-scales of high-z structures, that influence the 
initial evolution of the first objects, accordingly to the /nl values 
considered. On the other hand, it seems that feedback mechanisms 
are the key ingredient which largely wash out such original dis- 
crepancies and hinder the possibility of distinguishing the differ- 
ent models at later times (e.g. at z < 9). They, indeed, affect the 
amount of available cold star forming medium, and regulate the star 
formation activity, independently from the initial conditions. That 
is the reason why all the trends converge and are not tightly linked 
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to the non-Gaussian parameters. 

From our discussion it emerges that studies and observations of 
the primordial Universe can be very useful tools to try to ad- 
dress non-Gaussianities, much more powerful than low-redshift in- 
vestigations. In fact, signatures from first star formation episodes 
could be potential probes of primordial non-Gaussianities during 
the first Gyr of t he Universe. F or example, high-redshift, long 
gamma-ray bursts fcucchiara et al. 2011) seem to be highly cor- 
related with popIII st ar formation, as recently demonstrated by e.g. 
ICampisi et al.l(l201lh . So, they are suitable probes of the conditions 
of the Universe at early times, when non-Gaussian effects could be 
still visible in the baryon evolutio n. Similarly, first stars and galax- 
ies are detected up to z ^ 8 — 10 (iBouwens et al.ll20l7 ) and future 
instruments (like JWST) will likely see more and more of these ob- 
jects, providing important information about the status of the Uni- 
verse in its early infancy and possibly about the cosmic matter dis- 
tribution. Primordial episodes of star formation and metal enrich- 
ment from high-redshift objects have a very steep evolution (see 
Fig. [6]), and, even if the first onsets differ by only ^ 0.5 for 
/nl=0 and /nl=100, the corresponding metal filling factors dif- 
fer by ^ 2 orders of magnitudes. Despite the current difficulty of 
observing metal filling factors at such high redshift, they might pro- 
vide an interesting opportunity for future observational missions. 
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